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Abstract 

Vibrational spectra of long-range molecules are determined accurately and to arbitrary accuracy 
with the Canonical Function Method. The energy levels of the 0~ and 1^ electronic states of the 
^^Na2 molecule are determined from the Ground state up to the continuum limit. The method 
is validated by comparison with previous results obtained by Stwalley et al. [Q] using the same 
potential and Trost et al. [Q] whose work is based on the Lennard-Jones potential adapted to 
long-range molecules. 
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I. INTRODUCTION 



A new kind of high precision molecular spectroscopy is probing long-range forces between 
constituent atoms of molecules. This spectroscopy is based on using light to combine two 
colliding cold-trapped atoms into a tenuous molecule. 

The burgeoning field of "Photoassociation Spectroscopy" is allowing very precise measure- 
ment of lifetimes of the first excited states of Alkaline atoms and observation of retardation 
effects and long-range forces. It provides a means of probing accurately the weak interaction 
between these atoms 

The agreement between theory and experiment requires simultaneously a highly accu- 
rate representation of the interaction potential as well as a highly reliable method for the 
calculation of the corresponding energy levels. 

Since our aim is directed towards the latter problem, we make use of an alternative 
method to evaluate the energy levels for the potential at hand instead of comparing to the 
experimental values in order to assess the validity of our results. 

The determination of the vibrational spectra of these very tenuous molecules is extremely 
subtle specially for the highest levels which play an important role in photoassociation spec- 
troscopy. Thus a careful control of accuracy is needed in order to diagonalise the Hamiltonian 
without losing accuracy for all energies including those close to the dissociation limit. 

The magnitudes of potential energy, distance and mass values in these kinds of molecules 
stand several orders of magnitude above or below what is encountered in ordinary short- 
range molecules. 

For instance, the typical intramolecular potential well depth at the equilibrium distance 
of about 100 ao (Bohrs), is a fraction of a cm~^ while the reduced mass is several 10,000 
electron masses. All these extreme values require special numerical techniques in order to 
avoid roundoffs, divergences, numerical instability and ill-conditioning during processing. 

The method we use in this work adapts well to this extreme situation with the proviso 
of employing a series of isospectral scaling transformations we explain below. 

Since accuracy and its control are of paramount importance in this work, the canonical 
function method (CFM) is an excellent candidate because it bypasses the calculation of 
the eigenfunctions. This avoids losing accuracy associated with the numerical calculation 
specially with rapidly oscillating wave functions of highly excited states. 
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This method evaluates the full spectrum of the Hamiltonian, for any potential, to any 
desired accuracy up to the continuum limit. It has been tested succesfully in long-range and 
short-range potentials for atomic and molecular states. It describes faithfully bound and 
free states and is staightforward to code. 

This work is organised as follows: Section 2 is a description of the CFM and highlights 
the details of the method. Section 3 presents the results we obtain for the vibrational levels 
of the ^^Na2 molecule 0~ and 1„ electronic states. Section 4 is an additional validation of 
the method with the Lennard- Jones molecular potential used in a similar situation. We 
conclude in section 5. 

II. THE CANONICAL FUNCTION METHOD 

In this work we consider rotationless long-range diatomic molecules only. The associated 
Radial Schrodinger equation (RSE) is given by: 



where n is the reduced mass and V{r) is the potential energy of the interacting constituent 
atoms. 

The diagonalisation of the RSE is described mathematically as a singular boundary value 
problem. The prime advantage of the CFM is in turning it into a regular initial value 
problem. 

Primarily developed by Kobeissi and his coworkers 0] the CFM, consists, essentially, 
in writing the general solution as a function of the radial distance r in terms of two basis 
functions a{E] r) and l3{E] r) for some energy E. 

Picking an arbitrary point tq at which a well defined set of initial conditions are chosen 



§ i.e.: a{E;ro) = 1 with a'{E;ro) = and /3{E;ro) = with (3'{E;ro) = 1, the RSE is 



solved by progressing simultaneously towards the origin and oo. The effect of using different 
algorithms for the numerical integration will be examined in the next section. 

During the integration, the ratio of the functions is monitored until saturation signaling 
the stability of the eigenvalue spectrum 

One may define two associated energy functions: 




(1) 




a{E;r) 



(2) 
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and: 

The eigenvalue function is defined in terms of: 

F{E) = l^{E)-l^{E) (4) 

The saturation of the a{E] r)/f3{E; r) ratio as r progresses towards or oo yields a position 
independant eigenvalue function F{E). Its zeroes yield the spectrum of the RSE. 

Generally, one avoids using the wavefunction but if one needs it, the functions a{E; r) 
and l3{E] r) are used to determine the wavefunction at any energy with the expression: 

^{E; r) = ^{E; ro)a{E- r) + ^'{E; r^)(3{E- r) (5) 

where \l/(-E;ro) and '^'{E^tq) are the wavefunction and its derivative at the start point tq. 
The eigenfunctions are obtained for E = E^ where Ej. is any zero of F{E). 

The graph of the eigenvalue function has a typical tan(|i?|) shape versus the logarithm 
of the absolute value of the energy E as displayed in Fig. |l|. 
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FIG. 1: Behavior of the eigenvalue function F[E) with energy on a semi-log scale for the Vibrational 
spectra of the 0~ electronic state of the ■^^Na2 molecule. The vertical lines indicate the eigenvalue 
position. Energies are in cm^^. 
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III. VIBRATIONAL STATES OF THE ^^Na2 MOLECULE 



We apply the CFM to the calculation of the vibrational energy levels of a diatomic 
molecule where the interaction between the atoms is given by the Movre and Pichler potential 



We start with the 0^ electronic state of the ^^Na2 molecule. The corresponding potential 
is given by: 



V{r) = 2 [(1 - 3X) + Vl - 6X + 81X2] (6) 

where X = C(0g)/9r'^A. r is the internuclear distance and the parameter C'(Og) is such 
that: 

lim V{r) ^ (7) 

Identification of the large r limit yields the result: C3(0^) = C(0~)/3. We have used in the 
calculations below 6*3(0^) = 6.390 Hartrees.ap like |[T|]. The parameter A = 1.56512.10"^ Ry- 
dbergs is the atomic spin-orbit splitting. Given C^{Q^) and A, the equilibrium internuclear 
distance is = 71.6ao. 

We scale all energies with a factor Eq (usually cm"^) with the use of equation (1). Then 
we scale all distances with a typical length Lq transforming the RSE appropriately. This 
double transformation is reflected generally in the potential coefficients preserving thus the 
functional form of the potential. 

In order to gauge the accuracy of the spectra, we perform the integration of the RSE 
with two different methods: A fixed step Fourth-order Runge-Kutta method (RK4) and a 
Variable Step Controlled Accuracy (VSCA) method p. 

The VSCA method is based on a series expansion of the potential and the corresponding 
solution to an order such that a required tolerance criterion is met. Ideally, the series coef- 
ficients are determined analytically to any order, otherwise loss of accuracy occurs leading 
quickly to numerical uncertainties as discussed later. 

Table | shows the results we obtain with the RK4 method. The limitation of the RK4 
method to fourth order hampers the finding of levels beyond the 33rd (see Table 1). In order 
to find the higher levels we have to select an algorithm that enables us to tune the accuracy 
well beyond the fourth order. 

Pushing the accuracy within the framework of a fixed step method has the effect of 
reducing substantially the integration step. In order to avoid this problem, we use a variable 



TABLE I: Vibrational levels for the 0~ electronic state of the Na2 molecule as obtained with a 
fixed step RK4 method, Stwalley et al. ||T| results and the corresponding ratio. Levels 34-40 were 
not found by the RK4 method due to the precision limited to fourth order. 
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step that adjusts itself to the desired accuracy, the VSCA method. 

This method is powerful and flexible enough to find all the desired energy levels and 
allows us to find one additional level that was not detected before. It should be noted that 
the last three levels given by Stwalley at al. were extrapolated and not calculated. The 
agreement between our calculated levels and those of Stwalley et al. [Q] is quite good. We 
believe that the small discrepancy, increasing as we progress towards the dissociation limit, 
is due to a loss of accuracy associated with traditional methods in sharp contrast with the 
CFM. 

The estimation of accuracy of the results hinges basically on two operations, integration 
and determination of the zeroes of the eigenvalue function F{E). The superiority of the 
VSCA method is observed in the determination of the upper levels that are not detected 
by the RK4 method (see Tables 1 and 2). In addition, it is observed in the behavior of the 
eigenvalue ratio versus the index. While in both cases (RK4 and VSCA) the ratio increases 
steadily as the index increases because we are probing higher excited states, in the RK4 case 
it rather blows up as dissociation is approached. We use typically series expansion to order 
12 in VSCA with a tolerance of 10~^. In the root search of F{E), the tolerance required 
for a zero to be considered as an eigenvalue is 10^^^. This does not imply that we disagree 
as strongly as 0.13%, for instance, with the Ground state value (see Tables 1 and 2) found 
by Stwalley et al. [|I|] for the simple reason, we use a splitting energy A = 1.56512.10"^ 
Rydbergs corresponding to an equilibrium internuclear distance = 71.6ao. Stwalley et al. 
do not provide explicitly the value of A they use, and more recently Jones et al. P] provide 
a value that is slightly different. 

We treat next the 1^ electronic state of the ^^Na2 molecule. This state is higher that 
the 0~ and the number of vibrational levels is smaller because the potential is shallower as 
displayed in Fig. |. 

The potential associated with the 1„ electronic state of the ^^Na2 molecule is a lot more 
involved. Its VSCA implementation is particularly difficult because of the complex functional 
of the potential as we explain below. Analytically, the VSCA algorithm requires performing 
a Taylor series expansion to any order around an arbitrary point [Q. 

This is still, numerically, an open problem for arbitrary functions and the use of LISP 
based symbolic manipulation techniques produces quickly cumbersome expressions. Spe- 
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TABLE II: Vibrational levels for the 0~ electronic state of the Na2 molecule as obtained with 
Stwalley et al. results, the variable step controlled accuracy method (VSCA) method and the 
corresponding ratio. Levels 38, 39 and 40 of Stwalley et al. are extrapolated with LeRoy and 
Bernstein semi-classical formulae. 
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FIG. 2: Potential energy in cm~^ for the 0~ and !„ electronic states of the ^^Na2 molecule. The 
radial distance r is in ag units. 
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cial methods based on analytical fitting expressions are needed in order to turn the series 
coefficients into a more manageable form. 

The first step is to determine the 1„ electronic state of the ^^Na2 molecule by solving the 
Movre et al. secular equation such that: 

V{r) = A[-2^cos{^-^) - ^ - 1] (8) 

where a=-2-6X and X = C(lu)/9r^A. In addition, 6 = cos-\ }+^™^' ), and Q = i±f^. 
The parameter C{lu) is such that: 

lim V{r)^-^!^^ (9) 

Identification of the large r limit yields to the result: Cs{lu) = C(l„)(-\/7 — 2)/9. We have 
used in the calculations below C3(1m)=1.383 Hartrees.ag like Stwalley et al. The parameter 
A = 1.56512.10"^ Rydbergs is the same as for the 0~ state. 

Table |ITl| displays the results we obtain with the RK4 method that cannot find more than 
11 levels due to accuracy limitations. 

The next results for the 1„ electronic state of the ^^Na2 molecule are obtained with the 
VSCA method as shown in table ^V\. We find an additional 15-th level in contrast to Stwalley 
et al. who found fourteen and extrapolated the last two levels. 

The corresponding graph of the eigenvalue function is displayed in Fig. ^ below. 



IV. LENNARD-JONES MOLECULES 

We apply our methodology to the Lennard- Jones case. Our results are compared to the 
results obtained by Trost et al. . We start with the levels obtained with the RK4 method. 
The energy unit is e the depth of the potential well of the Asymmetric Lennard- Jones 
potential (ALJ): 

V{r) = C^{^)^-C2{lr (10) 
The Asymmetric-Lennard- Jones (ALJ) depends on Ci and C2 yielding an equilibrium 
distance at r = rmin and a potential depth of — e. Trost et al. [0] use the general parame- 
terisation: 
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TABLE III: Vibrational levels for the 1„ electronic state of the Na2 molecule as obtained with 
the RK4 method, Stwalley et al. results and the corresponding ratio. RK4 found only 11 levels 
and levels 15 and 16 of Stwalley et al. are found by extrapolation. 
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TABLE IV: Vibrational levels for the electronic state of the Na2 molecule as obtained by 

Stwalley et al., the variable step controlled accuracy method (VSCA) and the corresponding ratio. 

A new 15th level is obtained with the VSCA method. 

* Index Stwalley et al. (cm^^) VSCA (cm"^) Ratio 
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FIG. 3: Behavior of the eigenvalue function F{E) with energy on a semi-log scale for the 
electronic state of the ^^Na2 molecule. The vertical lines indicate the eigenvalue position. Energies 
are in cm~^. 
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It is scaled in such a way that the energy is expressed in units of the potential well depth 
e. When a = 6, /3 = 12 we obtain r^m = ^'^'^ ^ = ^c^- and the radial distance is 

in ^ where i? is a reduced scaled mass given hj B = Numerically Trost et al. 

use B = 10^ which is the order of magnitude encountered in long-range molecules within 
the framework of their system of units. The potential energy in these units is displayed in 
Fig. ^ The RK4 methods yields the results displayed in Table |V]. 
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FIG. 4: Potential energy of the Trost et al. model. The energy is in e units where e is the 
potential well depth and the radial distance is in where i? is a scaled mass. 

The accuracy limitation of the RK4 results in losing the uppermost level 24. Thus we 
move on to the results obtained with the superior VSCA method. All levels are obtained 
with the VSCA and the agreement with Trost et al. results is perfect as witnessed by the 
ratio values of Table 6. The eigenvalue function obtained with the VSCA method as a 
function of energy is displayed in Fig. |[ 



V. CONCLUSION 



The CFM is a very powerful method that enables one to find the Vibrational spectra 
of tenuous molecules where energies and distances are so remote from the ordinary short- 
range molecules case that special techniques should be developed in order to avoid numerical 
instabilities and uncertainties. 

The VSCA integration method used gives the right number of all the levels and the 
variation of the eigenvalue function F{E) definitely determines the total number of lev- 
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TABLE V: Quantum levels of a Lennard-Jones molecule in e units, the depth of the potential well 
as obtained by the RK4 method, Trost et al. and the corresponding ratio. Only the last 24th level 
was missed by the RK4 method. 
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TABLE VI: Quantum levels of a Lennard-Jones molecule in e units, the depth of the potential well 
as obtained by the VSCA method, Trost et al. and the corresponding ratio. 



* Index 


VSCA 


Trost et al. 




Ratio 


1 


-0.9410443 


-0.9410460 


1 


000002 


2 


-0.8299963 


-0.8300020 


1 


000007 


3 


-0.7276415 


-0.7276457 


1 


000006 


4 


-0.6336915 


-0.6336930 


1 


000002 


5 


-0.5478480 


-0.5478520 


1 


000007 


6 


-0.4698206 


-0.4698229 


1 


000005 


7 


-0.3992947 


-0.3992968 


1 


000005 


8 


-0.3359533 


-0.3359561 


1 


000008 


9 


-0.2794718 


-0.2794734 


1 


000005 


10 


-0.2295109 


-0.2295117 


1 


000003 


11 


-0.1857222 


-0.1857237 


1 


000008 


12 


-0.1477498 


-0.1477514 


1 


000010 


13 


-0.1152247 


-0.1152259 


1 


000010 


14 


-8.7766358(-02) 


-8.7766914(-02) 


1 


000006 


15 


-6.4982534(-02) 


-6.4982730(-02) 


1 


000003 


16 


-4.6469838(-02) 


-4.6469911 (-02) 


1 


000002 


17 


-3.1813146(-02) 


-3.1813309(-02) 


1 


000005 


18 


-2.0585953(-02) 


-2.0586161(-02) 


1 


000010 


19 


-1.2350173(-02) 


-1.2350373(-02) 


1 


000016 


20 


-6.6568735(-03) 


-6.6570240(-03) 


1 


000023 


21 


-3.0470500(-03) 


-3.0471360(-03) 


1 


000028 


22 


-1.0526883(-03) 


-1.0527480(-03) 


1 


000057 


23 


-1.9832170(-04) 


-1.9834000(-04) 


1 


000092 


24 


-2.6957891 (-06) 


-2.6970000(-06) 


1 


000449 


* 



16 




FIG. 5: Behavior of the eigenvalue function F[E) with energy on a semi-log scale for the Trost 
et al. Lennard-Jones molecule. The vertical lines indicate the eigenvalue position. Energies are 
in potential depth e units. 
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els. Generally it requires performing analytically Taylor series expansion to any order of 
an arbitrary potential function that might require the combination of numerical, symbolic 
manipulation and functional fitting techniques. Despite these difficulties, the results of the 
VSCA are rewarding. 

The use of the RK4 and the VSCA methods jointly paves the way to a precise comparative 
evaluation of the accuracy of the spectra obtained. In practice, the RK4 method can be 
used in optimisation problems whereas the VSCA is more adapted to the direct evaluation 
of the spectra. 

In this work, we did not consider molecular rotation, nevertheless the CFM is adapted 
to solve accrately Ro- Vibrational problems as well as any RSE diagonalisation problem. 

The methodology we describe in this work enables one to tackle precisely weakly bound 
states that are of great importance in low-energy scattering of atoms and molecules and 
generally in Bose- Einstein condensation problems 0. 
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